Luttinger sum rule for finite systems of correlated electrons 
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The validity of the Luttinger sum rule is considered for finite systems of interacting electrons, where the Fermi 
volume is determined by location of zeroes of Green's function. It is shown that the sum rule in the paramagnetic 
state is evidently violated within the planar t-J model at low doping while for the related Hubbard model, even 
in the presence of next-nearest-neighbor hopping, no clearcut exception is found. 
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I. INTRODUCTION 

In last decades extensive experimental and theoretical in- 
vestigations of high-temperature superconductors and other 
novel materials with strongly correlated electrons raised sev- 
eral fundamental questions regarding the Fermi liquid (FL) 
concept. One of the central assumptions in the FL theory is 
the existence of well defined Fermi surface with the volume, 
independent of electron-electron interactions. The latter prop- 
erty has firm basis in the Luttinger sum rule (LSR) 1 ' 2,3 , which 
represents the relation between the electron density and the 
poles or zeroes of the single-particle Green's function (GF) 
G(k, lo) at the chemical potential. 

To settle this issue several authors reconsidered the LSR, 
its validity and extensions regarding different aspects in con- 
nection with strongly interacting electrons. It has been shown 
that in one dimension (ID) LSR remains robust, although the 
usual FL is replaced by the Luttinger-liquid phenomenology 4 . 
LSR is valid even for models of strongly correlated electrons 
as the Kondo lattice model 5 . Recently, the LSR in the Mott 
insulating state has attracted attention. In this case the chem- 
ical potential is in the gap so there is no Fermi surface, how- 
ever the LSR can still apply for zeroes of GF and in this way 
the definition of Fermi surface is generalized to the concept 
of 'Luttinger' surface 6 and corresponding Luttinger volume 
(LV). In the insulator the LSR should remain valid for mod- 
els with the electron-hole symmetry as the Hubbard model 
with the nearest-neighbor hopping on a bipartite lattice 7 . On 
the other hand, extensions of the latter model can lead to 
the violation of the LSR in the insulating state 8 , or a novel 
form in ladder systems 9 . The LSR has to be also reformu- 
lated in the case of emergence of long range order 10 . Clearly, 
the most challenging question is the metallic state of strongly 
correlated electrons in the absence of any broken symmetry. 
There are indications from photoemission experiments that at 
least on some hole-doped cuprates the Luttinger volume is not 
conserved 11 ' 12 . Similar conclusions can be drawn from nu- 
merical studies of related relevant models as the t-J model 13 
and Hubbard model 14 near the Mott insulating state although 
deviations are modest. 

In this paper we consider the LSR on finite systems of cor- 
related electrons. We use the observation that the main ingre- 
dients of the theorem 1 are valid for finite systems, exploited 
rarely so far 15 . The underlying idea is that in this way one can 



directly test nontrivial models in different parameter regimes 
on the validity of LSR and the topology of the related Fermi 
surface and LV. While in this way it is still hard to prove the 
Luttinger theorem in the most relevant thermodynamic limit, 
it is much easier to show on its breakdown. As shown later, 
one should be careful in the classification of scenarios of LSR 
violation, since only some of them could remain relevant up 
to the thermodynamic limit. 



We investigate as examples two prototype models of 
strongly correlated electrons, most frequently studied in con- 
nection with superconducting cuprates, namely the two- 
dimensional (2D) Hubbard model and the t-J model on a 
square lattice. For comparison, we perform tests also for cor- 
responding ID models whereby the ID Hubbard model is ex- 
actly solvabl e) 16 ' 17 and the LSR should remain valid in the ther- 
modynamic limit. Within our numerical limitations presented 
results show that also 2D Hubbard model generally fulfills 
the LSR in the whole parameter range of the paramagnetic 
state, even in the presence of the next-nearest-neighbor hop- 
ping t' . Still, we find some exception which we classify as not 
clearcut. On the other hand, the 2D (and even ID) t-J model 
reveals evident deviations from LSR close to half filling. 



The paper is organized as follows. In Sec. II we briefly 
summarize the formalism underlying the Luttinger sum rule, 
in particular its relevance and application for finite-size sys- 
tems. We introduce the tight-binding models for interact- 
ing electrons, whereby the Hubbard model is the simplest 
one allowing the study of continuous development from the 
non-interacting case satisfying trivially LSR into the strong- 
correlation regime. The approach of avoiding the degeneracy 
with twisted boundary condition is described. We also clasify 
several possible scenarios of the LSR breakdown in finite sys- 
tems. In Sec. Ill we present numerical tests of the LSR for 
the Hubbard model and the t-J model, both on ID chains and 
2D square lattice. In particular, we try to pinpoint clearcut 
cases of the LSR breakdown, which are evident for 2D t-J 
model near half-filling. Sec. IV is devoted to conclusions and 
discussion. 
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II. LUTTINGER SUM RULE 
A. Formalism for finite systems 

We consider here the homogeneous system of interacting 
electrons on a lattice with periodic boundary conditions. In 
this case the LSRi^ can be expressed as 



N 



E 



i. 



(i) 



ks,G s (k,0)>0 



relating the number of fermions N to the value of zero- 
temperature T = single-particle GF G s (k, u — 0) at the 
chemical potential [i. General Tr 1 is in homogeneous case 
replaced by the sum over k and spin s. 

Let us recall few essential steps in the derivation of 
Eq. (Q} 1 ' 3 in order to see possible limitations and its proper 
application for finite systems. Generally one can express 



N 



(2) 



where <? s (k, wj) is the T > propagator at Matsubara uji — 
2tt(21 + I) /ft and = l/T (we use units h = k B = 1). 
With the definition of the self energy S s (k, Q — C — £ s(k) — 
1/G s {k, C) Eq. © can be rewritten as N = I x + I^, 



/ 2 = ^ L E / e s (k,o^s 5 (k,o« 



l 



where integration path T clockwise encloses real axis. The 
central point of the proof of the LSR is the observation that 
I2 in Eq. (0]) vanishes after limiting T — > 0. The latter is ar- 
gued with the construction of the functional Y', which is the 
contribution of all closed linked skeleton diagrams to thermo- 
dynamic potential. I2 can be represented as a full derivative 
of y', hence I2 = BY' /dde = 0^. We note that the existence 
and convergence of Y' is shown within the perturbation the- 
ory, which applies also for finite systems. It should be pointed 
out that deviations from LSR in concrete case discussed later 
on can be traced back to the nonvanishing I2 and in this way 
to the breakdown of perturbation theory for Y', 

In regular cases we have N = Ii which for T — > reduces 

to 



N 



^E 2Im [ £jC{^HG s (k,0)}]. (4) 



In finite systems we can express the GF for N electrons ex- 
plicitly in terms of eigenstates of systems with N— 1 . N, N+ 1 
particles, 



G s (k,C) 



E 

m 

E 
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(5) 



/ijv in Eq. (O is according to derivation 1 defined by the grand- 
canonical relation N = T9(lnf2) / 'd/j, in the limit T — > 0. 



The latter gives uniquely /in 



« +1 



E^ 1 )/2 under 



the provision that Eq is a concave function of N or at least 
E^ +1 + E^ 1 - 2E$ > 0, which is equivalent to the phase 
stability of the ground state (g.s.). It should be noted that in fi- 
nite systems the latter condition is not fulfilled in certain case 
which can be signature of a physical instability but as well 
just a finite-size effect. Another source of ambiguities in fi- 
nite size systems can arise from the degeneracy of the ground 
state I Oat), which can be partly removed by introduction of 
appropriate fields and twisted boundary conditions. 

It is easy to see that the contributions to Eq. (01 come from 
poles and zeros of G s (k, u> < 0). The latter appear in pairs, 
and give residuums ±1, respectively, mostly cancelling each 
other. As a rule we establish that finally the LSR, Eq. (fl}, 
is determined by the unpaired zero of GF at u> ~ which 
can appear either for u> > or u> < 0. This is a distinction 
to a macroscopic systems where the Fermi surface is located 
by the poles of GF at oj = . On the other hand, in a finite 
system the poles of GF are at ui ^ 0, since in Eq. (O in general 



- E» - E 



< (In and Eq +1 — Eq > /xjv due the choice 
of (In and concavity of Eq . Clearly, possible exceptions are 
when the latter condition is not satisfied as well when the g.s. 
is degenerate. 



B. Tight binding models 

In the following we test the validity of LSR for single-band 
tight-binding models for correlated electrons, 



H 



(6) 



where the prototype model is the Hubbard model with local 
repulsion 



Hint = U n 4i 



(7) 



and the second model is the t-J model discussed lateron. We 
consider finite systems with No sites on ID chain and 2D 
square lattice. Besides the nearest-neighbor hopping tij = t 
we investigate also the next-nearest-neighbor hopping ty = 
t', since the latter breaks the electron-hole symmetry, e.g. in 
the Hubbard model, which may affect the LSR 7 . The Hub- 
bard model has the advantage that by increasing U one can 
study continuous development from the reference system of 
noninteracting electrons (NIE) on a lattice, where the LSR is 
trivially satisfied, into a strong correlation regime where the 
breakdown of LSR can appear in several ways. It is evident 
that the interaction U can be treated within the standard per- 
turbation theory and therefore the arguments underlying the 
proof of LSR 1 should apply. 

We note that in finite systems even for U — a degeneracy 
of the g.s. can appear in the case of periodic boundary condi- 
tions due to discrete single -particle k = k,-,i = I, No. This 
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can be removed by introducing the twisted boundary condi- 
tions which are achieved by replacing — > exp^r.^ ) 
while the interaction in Eq. © is supposed not to depend on 
8. In this case, we are dealing on 2D square lattice with the 
reference system of NIE with the single-particle dispersion 

e(k) = — 2i(cos k x + cos k y ) — At 1 cos k x cos k y , (8) 

where k = k+8. To remove degeneracies, it is enough to take 
8 = (9 X , 6y) infmitesimally small, although also finite 8 may 
have a meaning, e.g. by minimizing the g.s. energy Eq(6). 



cases with even N,Nq. For smaller sizes, A^o < 8 one can 
evaluate G s (k, u>) via Eq. © finding all eigenstates using the 
exact (full) diagonalization within the basis states for given 
N, S z . To reach larger sizes, i.e. Ao = 16 for Hubbard model 
and A*o = 20 for the t-J model, respectively, we calculate 
the g.s. performing the exact diagonalization with the Lanc- 
zos technique, which is then also applied to calculate T = 
G s (k, u>) in a usual way 18 . It can be shown that such a proce- 
dure yields acute values for quantities of interest, in particular 
E^,^ N mdG 8 (k,0). 



C. Scenarios of sum-rule violation 

Before presenting results of our analysis, let us classify pos- 
sible ways of the LSR violation. We note that generally eigen- 
states of Eq.© can be sorted with respect to several quantum 
numbers: number of electrons N, total spin projection S z (as 
well as S), total momentum K, etc. According to the latter 
several scenarios are possible, e.g., for the Hubbard model, 
bearing in mind that the reference NIE system is paramag- 
netic: 

I) Turning on U > the character of the g.s. |0at) can change 
due to the crossing of levels with different symmetries. This is 
manifested by an abrupt jump of G(k, 0) and possible viola- 
tion of LSR. Such a case can be a signature of a macroscopic 
phase change, e.g., for S ^> 0^2 indicating a ferromagnetic 
instability (e.g., for U — > oo and A" = Ao ± 1 well known Na- 
gaoka instability), or merely a finite size effect (mostly change 
ofK). 

II) Similar, but more subtle, could be the effect of the level 
crossing in |0at + i) or |0at_i). In this case, /ijv is determined 
by the g.s. with AT + l,iV — 1 electrons, while G(k, 0), 
Eq. ((5]) does not necessarily contain matrix elements between 
N - 1, N, N + 1 g.s., if e.g. |5jv+i ~S N \> 1/2. As in sce- 
nario I such a level crossing could be a sign of a macroscopic 
instability. 

III) The g.s. of NIE can be degenerate on finite lattices. The 
latter can be eliminated by introducing small 8, but is removed 
also by U > 0. Such a case with possible level crossing be- 
tween different K at small U* (8) could lead to ambiguity, but 
is not the problem if the breakdown of LSR appears at larger 
U c > U*. 

IV) To avoid degeneracies of NIE, most rewarding are con- 
figurations with closed shells of electrons where we can fix 
8 = 0. Then most clearcut (excluding scenarios I, II) break- 
down of the LSR, Eq. (|T), can appear when by increasing 
U > U c one or several zeroes of G(k, ojq) cross the chem- 
ical potential, i.e. uj changes sign. 



III. NUMERICAL RESULTS 

In the following we present numerical results obtained for 
the ID and 2D Hubbard model and t-J model. In 2D lat- 
tices are chosen of the Pythagorean form No = L 2 + M 2 . 
To avoid complications in the interpretation we consider only 



A. Hubbard model 

ID systems: For the start we test the basic t' = Hubbard 
model, Eq. ©, on a chain with up to Ao = 14 sites. As ex- 
pected we find no violation of the LSR in cases with even N in 
the investigated range of < U < 30 t. Adding t' ^ breaks 
the particle-hole symmetry 7 and even at half-filling LSR could 
become questionable 8 . Nevertheless, choosing t' = —0.2 t 
we also do not find any evident breakdown following scenario 
IV of LSR at or away from half-filling. 

Several regimes with the violation of LSR (e.g., A^ = 6 
electrons on Ao = 14 sites for U > 6 t) can be mainly at- 
tributed to level crossing of |0jy+i) or |0jv— i) (scenario II). 
Also for half filled A^o = 12 system, where the degeneracy of 
the g.s. is removed by choosing small 8 we find N' = 10 
particles in LV instead of N = 12 above U c (8) < 4 t, 
which seems to represent the violation of type III. In this case 
N — 1, N, N + 1 g.s. preserve symmetry of NIE, however 
| Oat) has nonzero K. 

2D systems: t 1 = Hubbard model we investigate on square 
lattices with Ao = 8, 10, 16 sites. As a general rule (with 
some exceptions elaborated below), we find that the LSR re- 
mains satisfied for all considered N < No and moreover the 
topology of the LV remains that of NIE. In particular, we find 
no violation for half filling A^ = A^o, consistent with particle- 
hole symmetry^. 

Away from N = Nq one exceptional deviating case is 
N = 6 on Ao = 8 sites where at larger U > U c the LSR 
yields N' = 8, 10, 14 > N. Again, in this case NIE g.s. is 
degenerate and U c is 8 dependent and also for U < U c level 
crossing in A^ — 1, A^ + 1 g.s. is observed. All this indicates 
on combination of scenarios II and III. 

In our study of 2D t-t'-U model we use systems with 
Ao = 8, 10, 16 sites and t' = — 0.3t as frequently invoked 
for superconducting cuprates. We find no clear (type IV) vi- 
olation of LSR in range < U < 40 t and different fillings 
N < No- However, some deviations from LSR connected 
with level crossing scenario II are found at and away form 
half-filling even for closed shell configurations. E.g., A^ = 14 
electrons on A^ = 16 sites at U > 30 t result in N' = 20 
within the LV. Although |0jy) g.s. remains throughout S = 
and K = 0, in this case A^ + 1 g.s. changes from NIE value 
S = 1/2 to higher spin 5" > 1/2 at U > 30 t. 
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B. t-J model 

Another prototype model for strongly correlated electrons 
is the t-J model, 



a) NIE 



b) 0.1 t < J < 0.3 t 



c) s Cis + J ( S * ' S J _ \ n i n i)^ 

<ij> 



(9) 



where Ci S — Cj S (l — rii t - s ) are projected operators not al- 
lowing for a double occupancy at each site, and Sj are the 
local spin operators coupled with the exchange interaction J. 
The t-J model can be considered as the truncated Hubbard 
model at large U 3> t and J = At 2 /U. It is assumed that 
low-energy physics of both models is similar, therefore it is 
of interest whether they behave similarly with respect to LSR. 
We note that formally t-J model can be also written as the 
tight-binding model, Eq. (O, with the interaction 



H mt = J ^ (Si ■ Sj - jriirij) + 



(10) 



<tj> 



whereby we let U — » oo to get a projected model, Eq. (O. 
Such a formulation helps to formally connect the t-J model 
with finite J > continuously with NIE. Although the model 
is not perturbative in any limit, one can first consider the J > 
term as the perturbation and then performing the limit U — > 
oo. In this way it makes sense to compare the LSR within the 
t-J model directly to NIE with the dispersion, Eq.®. 
ID systems: Within the t-J model it makes sense to discuss 
cases N < No, i.e. representing the Mott insulator doped with 
Nh = No — N holes. The model on a chain at J — ► behaves 
as the U = oo Hubbard model obeying LSR as confirmed 
in our numerical calculations. Finite J > is nontrivial and 
some breakdowns of LSR are found at low doping Nh <C Nq. 
For No — 16 and Nh = 2, N = 14 representing a closed shell 
configuration we find at J > 0.7 t N' — 24. Since the |0jv) 
has S = 0,K = Oas well as N - 1,N + 1 g.s. have S = 1/2, 
this case represents a clean violation of LSR of type IV. 
2D systems: Here, one should also monitor the concavity of 
Epj which can be violated at low doping, but not for even 
TV considered in our study. In general, results and conclu- 
sions for the LV within the t-J model are quite similar to the 
corresponding large- U Hubbard model. Essential differences 
appear at higher parameter J > 0.1 1 and at Nh <C Nq. We 
investigate further on lattices with Ao = 16, 18, 20 sites and 
Nh < 6. 

The most clear counterexample of LSR following scenario 
IV is the closed-shell system with Nh = 2 on No = 20 sites. 
The corresponding LV for NIE is shown Fig. la, together with 
wave-vectors appearing onaJVo = 20 lattice. On the other 
hand, we find that for 0.1 t < J < 0.3 t the LV contains 
A' = 22 particles, as shown in Fig. lb. Here, |0jv) nas S = 
0, K = 0, and no level crossing in N — 1, N + 1 g.s. with 
S = 1/2 is detected in this regime, excluding scenario II. For 
larger J/t even more drastic violations are found. For J = 
0.4 t,J= 0.5 t we find N' = 24 and N' = 32, respectively, 
as represented in Figs.lc,d. It is expected that introduction of 
t' should even increase the deviation from the LSR 19 . 



n,n ) 



c) J = 0.4 t 




7r,jr) 



J0,'<k 



\ . /o 
o\ / • 



T,0) 



Figure 1 : Luttinger volume for closed shell configuration of N — 18 
electrons on iVo = 20 square lattice for: a) noninteracting electrons, 
b) t-J model with 0.1 1 < J < 0.3 t, c) J = 0.4 t, and d) J = 0.5 t. 
Gray (blue) circles represents (nonequivalent) wave vectors k outside 
(inside) the Luttinger volume. 



Results for the violation of the LSR, presented in Figs, lb-d 
are not unexpected and qualitatively consistent with other (in- 
direct) indications of the violation of the LSR within the t-J 
model, obtained via the high-T expansion of momentum dis- 
tribution rtk 13 or a straightforward interpretation of spectral 
functions A(k, u>), calculated on small systems 19 . Namely, in 
the range of parameters for cuprates J ~ 0.3 t results within 
the t-J model show too large LV, N' > N, whereby the vio- 
lation is modest 13,19 , the conclusion is consistent with experi- 
mental finding for hole doped cuprates l M2k. 

In a doped system larger J plausibly enhances the antifer- 
romagnetic order (remaining always finite-range in small sys- 
tems) which can lead to the formation of hole pockets, which 
should show up in small Fermi surface and corresponding LV 
forming around k = (tt/2, 7t/2). Indeed such tendency is ev- 
ident in Figs. lc,d where at largest J = 0.5 t only the point 
closest to k = (7r/2, 7t/2) remains outside the LV. 



IV. CONCLUSIONS 

In conclusion, we have shown that the investigation of 
finite-size systems can provide a non-trivial test of the LSR 
and its validity for strongly correlated systems. Several prop- 
erties of finite systems, in particular the degeneracy of g.s., 
can complicate the interpretation of results. Nevertheless, the 
cases with the major violation of the LSR should be easily de- 
tectable also in small systems. In this study we examine only 
prototype Hubbard and t-J models, still the application of the 
method to other models of correlated electrons is straightfor- 
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ward. 

It should be pointed out that in all evident (type IV) cases 
the breakdown of the LSR appears by a continuous variation 
of GF zero ujo ~ across the chemical potential, which in a fi- 
nite system is connected with a divergence of the correspond- 
ing S(k, uj ~ 0) — > oo. Such a transition remains an evident 
possibility in an insulator even in the limit N — > oc 8 . On the 
other hand, in a metal (where in a normal FL the Fermi surface 
is determined by poles of GF) such a scenario would represent 
a major modification of the LV concept and therefore the limit 
N — > oo allowing for several nontrivial scenarios should be 
considered with care. 

For the Hubbard model on ID chain and 2D square lattice, 
both for t' = and t' ^ 0, we do not find a clearcut viola- 
tion (type IV) of the LSR. While this is rather a test of our ap- 
proach for ID systems where exact results 16, 17 indicate that for 
N — > oo the LV remains that of the NIE with the singularity at 
kp = ttN/2Nq away from half filling. For 2D system the va- 
lidity of the LSR would represent the confirmation of the fact 
that the g.s. even at large U > is adiabatically connected to 
that of NIE, as far it represents the paramagnetic state. Sev- 
eral examples of type II, III deviations still point on a cautious 
interpretation, since systems studied numerically are small. 
In particular for 2D lattice, large separation of discrete levels 
for NIE in considered systems indicate that it would be hard 
to observe small deviations from LSR. Another possibility in 
D > 1 is a change of the Luttinger-surface form while pre- 
serving LV, allowed and expected for U > 0. In our approach 
the latter would show up with few zeroes of GF entering and 
leaving the LV. Since nonequivalent zeros can hardly cross fx 



at the same parameter in a finite system, this would lead to 
intermediate violation of LSR. Again, our systems seem to be 
too small to observe this phenomenon. 

Within the t-J model we observe the violation of the LSR, 
particularly evident for the 2D systems corresponding to low 
doping of the Mott insulator as relevant to superconducting 
cuprates. This can be interpreted the g.s. of the model not 
being adiabatically connected to the g.s. of NIE, although for- 
mally the latter would be allowed by performing the perturba- 
tion in J > and letting U — > oo in Eq. dlOb . It is not sur- 
prising that the latter extrapolation can break the connection 
to NIE. In our study we follow numerically only J > 0.1 t. 
On the other hand, it is clear that J — > case should be equiv- 
alent to the Hubbard model with U — » oo. In the latter regime, 
the Hubbard model (or t-J model with J — > 0) close to half- 
filling can show instability towards partly or fully spin po- 
larized states S ^> as well as more pronounced finite-size 
effects, which are only partly studied in the present work. 

Our results are consistent with other theoretical indications 
for the deviations from LSR 13 '— as well as experiments on 
cuprate o 11 ' 12 . Still the speculations on the origin and the ex- 
trapolation to a macroscopic system are delicate. While the 
deviation from the LSR could indicate a general breakdown 
of the Fermi liquid concept or be a sign of inherent insuffi- 
ciency of the model, the deviations could as well disappear in 
the limit N — > oo in a metallic system while persisting (or 
not) e.g. in a Mott insulating state. 

One of the authors (P.P.) thanks M. Potthof for useful sug- 
gestions regarding the LSR. 
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